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ABSTRACT 

A hierarchical Bayesian method is applied to the analysis of Type-IA supernovae (SNIa) ob- 
servations to constrain the properties of the dark matter haloes of galaxies along the SNIa 
lines-of-sight via their gravitational lensing effect. The full joint posterior distribution of the 
dark matter halo parameters is explored using the nested sampling algorithm MultiNest, 
which also efficiently calculates the Bayesian evidence, thereby facilitating robust model 
comparison. We first demonstrate the capabilities of the method by applying it to realistic 
simulated SNIa data, based on the real 3-year data release from the Supernova Legacy Survey 
(SNLS3). Assuming typical values for the halo parameters in our simulations, we find that a 
catalogue analogous to the existing SNLS3 data set is incapable of detecting the lensing signal, 
but a catalogue containing approximately three times as many SNIa does produce robust and 
accurate parameter constraints and model selection results for two halo models: a truncated 
singular isothermal sphere (SIS) and a Navarro-Frenk-White (NFW) profile, thereby validat- 
ing our analysis methodology. In the analysis of the real SNLS3 data, contrary to previous 
studies, we obtain only a very marginal detection of a lensing signal and weak constraints on 
the halo parameters for the truncated SIS model, although these constraints are tighter than 
those obtained from the equivalent simulated SNIa data set. This difference is driven by a 
preferred value of rj w 1 in the assumed scaling-law a oc If between velocity dispersion 
and luminosity, which is somewhat higher than the canonical values of r} = j and rj = | for 
early and late-type galaxies, respectively, and leads to a stronger lensing effect by the halo. 
No detection of a lensing signal is made for the NFW model. 
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1 INTRODUCTION 

In using Type la supernovae (SNIa) as 'standardizable' can- 
dles to constrain cosmological parameters, one typically assumes 
that the universe is homogeneous and isotropic, and therefore 
one ignores gravitational lensing effects due to cosmic struc- 
ture along the line-of-sight to each SNIa. In conventional cos- 
mological SNIa analyses, this effect is usually regarded as an 
additional source of uncertainty, which adds ex tra scatter to the 
brightness of SNIa that increases with redshift dKantowski et al.l 
19951 lFriemadll99dlWambsganss et al.lll997llHolz & Walj|l998l, 
Bergstrom et alj|2000h . Fortunately, owing to flux conservation, 



the effects of gravitational magnification and demagnification av- 
erage out and are therefore expected to lead to negligible bias 
in cosmological pa rameter estimates (see, e.g., ISarkar et alj|200a . 
IJonsson et al]|2008h . 
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Nonetheless, the gravit ational lensing of SNIa can itself be 
used to constrain cosmology jMetcalJl9 99. Dodel son & Vallinottol 
2006, Zentner & Bhattacharya 2009) and/or the pr operties of the 
lensing matter dRaucdll99lUMetcalf & Silklll999h . In the latter 
case, one performs a complementary analysis to cosmological pa- 
rameter estimation by instead assuming a particular background 
cosmological model (i.e. fixing the cosmological parameters to 
some concordance values) and using the observed distance moduli 
to constrain the nature of the cosmic structure, such as the proper- 
ties of dark matter haloes, along the lines-of-sight to the SNIa. In 
principle, one might even hope to perform a joint analysis to con- 
strain the background cosmological parameters and the nature of 
the cosmic structure simultaneously, but such an approach is likely 
to suffer from strong degeneracies between parameters. 

An early t entative detect i on of gravitational lensing of SNIa 
was made by IJonsson et al. I d2007tj using a sampl e from the 
Great Observatories Origin s Deep Surv e y (GO ODS; iRiess et ail 
|2004 IStrolger et ail |2004 IRiess et al.l UOOlb . More recently, 
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Kron borg et all d2010h focussed on the detection of a gravitational 
lensing signal by assuming properties for dark matter haloes and 
repor ted a positive result at the 99 per cent confidence level. More- 
over, IjQnsson eTa l, (2010a) used 24 high-redshift (0.4 < z < 
1.8) SNIa from GOODS to constrain the properties of dark mat- 
ter haloes of ga laxies also contained within GOODS. This study 
was extended in ljonsson et alj ( feOlObl) by using 175 high-redshift 
(0.1 < z < 1) SNIa fro m the 3-year data r elease of the Supernova 
Legacy Survey (SNLS3; lAstier et al.ll200q) to constrain the haloes 
of galaxies in the deep Canada-France-Hawaii Telescope Legacy 
Survey (CFHTLS) fields. Although the SNIa in SNLS are typically 
not as distant as those from GOODS, they are far more numerous 
and selected in a more homogeneous way. Jonsson et al. (2010b) 
report the detection of a gravitational lensing signal at the 92 per 
cent confidence level. 

In this paper, we also use high-redshift SNIa from SNLS3 
to constrain the properties of dark matter haloes of galaxies in 
the CFHTLS fields that intersect the SNIa lines-of-sight. Our sta- 
tistical methodology differs greatly, however, from that used by 
Ijonsson et al. I<201PJ|B and ot her previous studies. As recently dis- 



cussed by Mar ch et al. U201 ]]), the usual x z -method used to con- 



strain cosmological parameters and/or the nature of c osmic struc- 
ture from lightcurve fits to SNIa observations ( see e.g. lAstier et al.1 
120061: iKowalski et aT]|2008l; IConlev et alJlioT 1) suffers from some 
several shortcomings in terms of its statistical foundations and ro- 
bustness, including not allowing for rigorous model checking and 
not providing a reliable framework for the evaluation of system- 
atic uncertainties. Consequently, we instead analyse the SALT- 
II lightcurve fits of the SNIa observations using the statistically- 
principled and rig orous Bayesian hierarchical method (BHM) of 
iMarch et al.l d201lh to obtain a robust effective likelihood function 
giving the probability of obtaining the observed SNIa data (i.e. the 
parameter values obtained in the S ALT-II lightcurve fits) as a func- 
tion of the parameters of the dark matter halo models assumed for 
the galaxies along the lines-of-sight to the SNIa. Moreover, rather 
than exploring the parameter space of t he dark matter halo m odels 
using simple gridding methods (see e.g. ljonsson et alj2010al lbh. we 
instead sample from the full joint posterior distribution of the dark 
matter halo parameters using a nested sa mpling algorithm (Skill ina 
l2004lFeroz & HobsonlfeoOSl i lFeroz et al.ll2009f) . This enables us to 
explore all the halo model parameters simultaneously and allows 
for straightforward marginalisation over subsets of them. The al- 
gorithm also efficiently calculates the Bayesian evidence, thereby 
facilitating robust model comparison. 

The outline of this paper is as follows. In Section 2 we give 
a brief summary of Bayesian inference methods, followed in Sec- 
tion 3 by a description of our Bayesian methodology for using grav- 
itational lensing of SNIa to constrain the properties of dark matter 
haloes of the galaxies intersecting the lines-of-sight. In Section 4, 
we describe the SNLS3 supernovae and galaxies data used in our 
analysis. We test the performance of our Bayesian methodology in 
Section 5 by applying to realistic simulated data based on the real 
SNLS3 data, before analysing the real data sets and presenting our 
results in Section 6. We give our conclusions in Section 7. 

Throughout the paper, we assume a spatially-flat concordance 
ACDM background cosmology, characterised by the parameters 
C = {fl m ,a, n A ,o,H } = {0.27, 0.73, 0.7}. Finally, we note that 
this pa per may be cons idered as complementary to our companion 
paper {March et al.l2012h , in which we use the BHM to analyse the 
SNLS3 catalogue, together with additional low-redshift SNIa, to 
constrain the background cosmological model, assuming no gravi- 
tational lensing along the lines-of-sight to the SNIa. 



2 BAYESIAN INFERENCE 

Our analysis methodology is built upon the principles of Bayesian 
inference, which provide a consistent approach to the estimation of 
a set of parameters © in a model (or hypothesis) H for the data D. 
Bayes' theorem states that 



Pr(0|D,H) 



Pr(D|0,#)Pr(0|#) 



(1) 



Pr(D|ff) 

where, for brevity, we denote Pr(0|D,i?) = P(&) as the pos- 
terior probability distribution of the parameters, Pr(D|0,i/) = 
£(0) as the likelihood, Pr(&\H) = 7r(0) as the prior, and 
Pr(D|ff) = Z as the Bayesian evidence. 

In parameter estimation, the normalising evidence factor is 
usually ignored, since it is independent of the parameters 0, and 
inferences are often obtained by taking samples from the (un- 
normalised) posterior using standard MCMC sampling methods, 
where at equilibrium the chain contains a set of samples from the 
parameter space distributed according to the posterior. This poste- 
rior constitutes the complete Bayesian inference of the parameter 
values, and can be marginalised over each parameter to obtain in- 
dividual parameter constraints. 

In contrast to parameter estimation problems, for model se- 
lection the evidence takes the central role and is simply the factor 
required to normalize the posterior over 0, 



(2) 



Z = / £(0)7r(0)d"0, 



where D is the dimensionality of the parameter space. As the av- 
erage of the likelihood over the prior, the evidence is larger for a 
model if more of its parameter space is likely and smaller for a 
model with large areas in its parameter space having low likeli- 
hood values, even if the likelihood function is very highly peaked. 
Thus, the evidence automatically implements Occam's razor. The 
question of model selection between two models Ho and Hi can 
then be decided by comparing their respective posterior probabili- 
ties given the observed data set D, as follows 



Pr(ij r ilD) _ Pr(D[ij'i)Pr(H 1 ) _ Z x Vx(H Y 



(3) 



Pr(ff |D) Pr(D|tf )Pr(Ho) Z Pr(H 

where Pr(Hi)/ Pt(Hq) is the a priori probability ratio for the two 
models, which can often be set to unity but occasionally requires 
further consideration. 

Evaluation of the multidimensional integral in Eq. l(2j is a 
challenging numerical task. Standard techniques like thermody- 
namic integration are extremely computationally intensive which 
makes evidence evaluation at least an order of magnitude more 
costly than parameter estimation. Some fast approximate meth- 
ods have been used for evidence evaluation, such as treating the 
posterior as a multivariate G aussian centered at its peak (see e.g. 
Hobson & McLachlanll2003h . but this approximation is clearly a 
poor one for multimodal posteriors (except perhaps if one performs 
a separate Gaussian approximation at each mode). The Savage- 
Dickey density ratio has also been proposed (see e.g. lTrottall2007]) 
as an exact, and potentially faster, means of evaluating evidences, 
but is restricted to the special case of nested hypotheses and a 
separable prior on the model parameters. Various alternative infor- 
mation crite ria for astrophysical model selection are discussed by 
iLiddld ( I2007I) . but the evidence remains the prefer red method. 

The nested sampling approach, introduced bv lSkillind J2004h . 
is a Monte Carlo method targeted at the efficient calculation of the 
evidence, but also prod uces posterior infere nces as a by-product. 
iFeroz & HobsorJ d2008l) and lFeroz et"al] d2009l) built on this nested 
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sampling framework and have introduced the MultiNest algo- 
rithm which is very efficient in sampling from posteriors that may 
contain multiple modes and/or large (curving) degeneracies and 
also calculates the evidence. This technique greatly reduces the 
computational cost of Bayesian parameter estimation and model 
selection and has alrea dy been applied to a number of problems 
in astrophysics (see e.g.lFeroz & Hobsonll201ll ; iFeroz et al]|201ll ; 
Bridg es et al]|2009l : IFeroz et al]|2010l) . We employ this technique 
in this paper. 



3 ANALYSIS METHODOLOGY 
3.1 Definition of the SNIa data 

In practice, there are no perfect astronomical standard candles. 
In particular, SNIa have absolute magnitudes that vary by about 
±0.8 mag in the B-band due to physical differences in how each 
supernova is triggered and also due to absorption by its host galaxy. 
Nonetheless, SNIa do constitute a set of 'standardizable' can- 
dles, since by applying small corrections to their absolute mag- 
nitudes, derived by fitting multi-wavelengths observations of their 
lightcurves, one can reduce the scatter considerably, to around 
±0.15 mag in the B-band. In essence, SNIa with broader light 
curves and slower decline rates are intrinsically b righter th an those 
with narrower light curves and fast decline rates ( Phillips 1993). 

Several methods are available for fitting SNIa lightcurves 
(and constraining cosmological parameters), including, amongst 
others, the M ulti-Colour Lightcurve Shape (MCLS) strateg y 
jjha et alj|2007t) , CMAGIC jWang et alJliooalConlev et alfeOOfl) . 
and the Spectrall y Adaptive Lightcurve Template (SALT) method 
dGuv et alj [2007). the current version of which is SALT-II. The 
relative merits of these methods is a topic of much debate, but the 
SALT-II method is particularly attractive for our purposes, since 
(unlike MCLS) it first fits each SNIa lightcurve to obtain three pa- 
rameters controlling the SN magnitude, stretch and colour correc- 
tions to a template 'learned' from nearby and distant SN; only in a 
second, separate step a re these fits used to constrain cosmological 
parameters. Following Ma rch et al. we may therefore use 

the products of the first step as the inputs to a statistically rigorous 
Bayesian hierarchical model. 

Our analysis takes place after the selection cuts, lightcurve fit- 
ting and Malmquist correction have been performed. For each se- 
lected SNIa, in addition to an estimate I of its redshift and an asso- 
ciated uncertainty cr z , derived from observations of its host galaxy, 
we take as our basic data the output from the SALT-II lightcurve 
fitting algorithm, which produces the best-fit values: m* B , the rest 
frame B-band apparent magnitude of the supernovae at maximumn 
luminosity; X\, the stretch parameter related to the width of the fit- 
ted light curve; and c, the colour excess in the B-band at maximum 
luminosity. These are supplemented by the covariance matrix of the 
uncertainties in the estimated lightcurve parameters, namely 



covariance matrix C%. The 'observed' distance modulus /i° 
each SN is then modelled as 



C = 



xi, c 
2 



Therefore, our basic input data for each SN (i = 1, . . . , Nsn) are 



(4) 



(5) 



and we assume (as is implicitly the case throughout the SNe lit- 
erature) that the vector of values {rh* B:i , Ci) for each SN is 
distributed as a multivariate Gaussian about the true values, with 



obs ~ * 

Mi = m s, 



Mi + ax lyi — /3£j, 



for 



(6) 



where Mi is the (unknown) B-band absolute magnitude of the SN, 
and a, f} are (unknown) nuisance parameters (assumed the same 
for all SN) controlling the stretch and colour corrections. 

It should be noted that a sophisti cated Bayesian hierarchica l 
method has recently been proposed by iMandel et"ai] J2009l LlQlch 
to fit optical and infrared lightcurve data. This may provide a more 
robust technique for defining the basic SNIa data that we use in our 
subsequent analysis, but we leave the investigation of this issue to 
a future work. 



3.2 Computing the predicted distance moduli 

In using SNIa to constrain cosmic structure, the predicted dis- 
tance modulus must include the effect of gravitational magnifica- 
tion due to cosmic structure along the line-of-sight to each super- 
nova. In particular, we assume this magnification is due to dark 
matter haloes associated with known galaxies intersecting the line- 
of-sight. Thus, for each SNIa, the predicted magnification depends 
on the sets of parameters {g, h}, where 



9 = {ZgabflgabMsjT 1 , 



"gal 



Eal 0^ l ,Mg sal ,r"^} (7) 



contains the redshift, sky position, absolute B-band magnitude and 
spectral type of the N sa i galaxies that intersect the line-of-sight to 
the SNIa, and h contains the parameters of the assumed dark mat- 
ter halo model for these galaxies (see Section [3~3l >, In general, the 
parameters g are (naturally) different for each SN and are assumed 
known; one wishes to place constraints on the unknown halo pa- 
rameters h, which (perhaps unrealistically) are assumed common 
to all the foreground galaxies. 

To compute the predicted magnific ation of a SNIa, we u se the 
weak-lensing approximation (see, e.g.. Schneider et al. 1992), the 
validity of which was checked by Ijonsson et al] feOlObl) by com- 
paring it with a ray-tracing algorithm. For the SNIa sample we 
consider here (see Section |4~Tt , which have relatively low redshifts, 
the weak-lensing approximation was found to be accurate to within 
~ 5 per cent. In this approximation, the predicted distance modu- 
lus of a SNIa, expressed in terms of magnitudes, is related to the 
convergence k along its line-of-sight by 

H(z,C,g,h) w no{z,C) - 2.17[« lo s(ff, fc) - «*>(&)], (8) 

where no(z,C) is the predicted distance modulus for our as- 
sumed cosmological parameters C, neglecting gravitational lensing, 
n\ os (g, h) is a sum over the contributions to the convergence from 

each galaxy along the line-of-sight, so that K\ as = YljS 1 K g a i> an d 
Ki,(h) represents the compensating effect of the background den- 
sity and acts as a normalisation allowing the magnification relativ e 
to a homogeneous universe to be computed (lon sson et alj|2010bh . 
Flux conservation implies that (k) = and to ensure this condition 
is satisfied we set Kb(h) = {KiosQ?) h)), where the latter is calcu- 
lated using a large number of randomly selected lines-of-sight. 



3.3 Halo models 

Our primary goal is to constrain the parameters h that describe the 
properties of the dark matter haloes associated with the galaxies 
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along the lines-of-sight to the SNIa. We assume that each galaxy 
contributes to the convergence by an amount 

s(€) 



(9) 



where the surface density, £(£), is obtained by projecting the mat- 
ter distribution onto a lens plane, 



(10) 



where £ is a vector in the plane and y is a coordinate along the line 
of sight; it is at this point that the density profile p(r) of the dark 
matter halo enters the calculation. The denominator in ® is the 
critical surface density, 

E c =^i, (11) 



4ttG A As ' 

which, in turn, depends on the angular diameter distances between 
the observer and the source, A, the observer and the lens, A, and 
the lens and the source, As- These distances are computed from 
the redshifts of the SN and galaxy respectively, assuming our con- 
cordance background cosmology. 

In this paper, we investigate two different models for the den- 
sity profile p(r) of the dark matter halo: t he truncated singula r 
isothermal sphere (SIS) and the NFW profile dNavarro et all 19971) . 
both of which we now describe. 



3.3. 1 Truncated singular isothermal sphere 

The radial density distribution of a singular isothermal sphere (SIS) 
is given by 



p(r) = 



(12) 



2-kG r 2 ' 

which depends on the single free parameter a, the one-dimensional 
velocity dispersion of its constituent particles. Since the SIS profile 
has a divergent total mass, we truncate it at a radius r = rt, which 
is thus a second free parameter. 

The surface density of a truncated SIS is easily calculated us- 
ing l llOt and found to be 







arctan \J r 2 /£ 2 — 1 if £ < rt 



(13) 



if € > rt, 

which can, in turn, be substituted into (|9} to obtain the convergence 
K gai (£) due to the galaxy dark matter halo. 

To allow for and investigate the relationship b etween galaxy 
lumino sity and velocity dispersion, we follow | jonsson etaU 
l l2010bl) and adop t the Faber- Jackson {Faber & JacksorJ 19761) and 
Tully-Fisher-like jTullv & Fishenl 19771) scaling law 

l y 

a) ' (14) 

where A is a fiducial luminosity, which we take to be 
A = 10 10 h~ 2 Lq in the B-band, and <r» is the corresponding 
fiducial velocity dispersion. In terms of absolute B-band magni- 
tudes, with which we will be working, the scaling relation becomes 



a = ff*10 



-1(«b-Mb)/ 2 -5 



(15) 



where M B = —19.52 + 5 log 10 h and h is the Hubble constant in 
units of 100 km s _1 Mpc -1 . The (aperture) mass-to-light ratio of 
the galaxies is determined by 77, since 

M(r < R) _ 2ol{L/L t ) 2rl R 



GL 



oc L 



->>)- 



(16) 



provided R ^ r t . Thus, for example, one has a constant mass-to- 
light ratio M/L oc L° if 77 = 0.5. 

Fo r the truncation radius, we again follow Ijonsson et al.l 
(2010b) and assume it obeys a scaling law of the form 



(17) 



where r» is a fiducial truncation radius. Since we only include 
galaxies located a distance 8 C from the position of the SN la, the 
truncation radius has an effect only if r t / Ad < C . For 9 C — 60" 
and Zg a i < 1 this corresponds to r t < 300/i -1 kpc. 

Thus, for the truncated SIS halo model, we wish to constrain 
the four parameters h — {7, 77, <r*, r,}. 



3.3.2 Navarro-Frenk-White profile 

The radial density distribution of the Navarro-Frenk-White (NFW) 
profile is given by 

Sc pc 



p{r) = 



(r/r s )(l + r/r s 



(18) 



where p c — 3H 2 (z) /8nG is the critical density for closure of the 
universe at the redshift, z, of the halo and H (z) is the Hubble pa- 
rameter at that same redshift. The scale radius r s = r2oa/c is a 
characteristic radius of the cluster, where c is a dimensionless num- 
ber known as the concentration parameter, and 

«. = gg ° 3 09) 



3 ln(l + c)-c(l + c)- 1 

is a characteristic overdensity for the halo. The virial radius, r2oo, 
is defined as the radius inside which the mass density of the halo is 
equal to 200p c - The profile therefore depends on two free parame- 
ters: the virial radius r^oo and the concentration parameter c. The 
mass of an NFW halo contained within the radius r2oo is 

m - <- \ 8007r 3 8007r P™ 3 r»m 
M200 = M(r ^ r 2 oo) = — — p c r 200 = — — r 20 o, (20) 

where p m is the mean matter density of the universe and Q m is the 
matter density parameter at the redshift z of the halo. 

The surface density of the NFW profile is given by 
teartelmaniill 19961) 



s(0 = 



2r s 5 cPc 

2r s <5 c p c 
3 

2r s 5 aP? 



. 2 arctanh •> / i-^- 
yfi^ V 1+x 



arctan . 



(X<1) 
(x=l) 

(x>l) 



(21) 

where x = £/r s is a dimensionless projected radial distance in the 
lens plane. Once again, this expression may be substituted into ^ 
to obtain the convergence /t ga l(£) due to the galaxy dark matter 
halo. 

By analogy with the Tully-Fisher and Faber-Jackson rela- 
tions, to allow for and investigate the relationship between the 
galaxy luminosity and the virial radi us in a simil a r man ner to 
the truncated SIS profile , we follow Ijonsson et al. (2010b) and 
Kleinheinr ich et all d2004l) and adopt the scaling law 



r-200 = J-200 



L 



(22) 



where, once again, A is a fiducial luminosity that we take to be 
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L* = 10 10 h 2 Lq in the B-band, and rJ o is the corresponding 
fiducial virial radius. The virial mass-to-light ratio is thus given by 

M 200 8 00tt rf oc 
= n ~ — ~L~ 



Nuisance parameter 



Symbol 



Value 



oc L 



37J-1 



(23) 



L 3 

For example, a constant mass-to-light ratio M/L oc L° requires 
A=±. 

Thus, for the NFW halo model, we wish to constrain the three 
parameters h = {A, c, r* 00 }. 

3.4 Likelihood function 

To construct the likelihood function for the SNIa data, w e adop t 
the Bayesian hierarchical model proposed bv lMarch et alJ j201lh . 
This has been shown to deliver tighter constraints on cosmological 
parameters than the usual y 2 -method, and simultaneously provides 
a robust statistical framework for the full propagation of systematic 
uncertainties to the final inferences. 

For the ith SN, we require the likelihood of the input data l[5j 
given the parameters of our model, namely 



Pr(m-s,i, St\,i, Cj, Zi\h, a, 0, a int ), 



(24) 



which also depends on the cosmological parameters C, the 
covariance matrix d of the uncertainties on the input data 
(m* B i ,xi t i, Ci), and the uncertainty a z .i in the estimated redshift 
Zi, all of which are assumed known. In particular, we seek to con- 
strain the unknown halo parameters h, global colour and stretch 
correction multipliers a and /3, and the intrinsic dispersion oi nt of 
SNIa absolute magnitudes (all of these are assumed to be 'global' 
parameters, i.e. common to every SN). Finally, the full likelihood 
function is give n by the product of the likelihoods J24t for each SN. 

Following lMarch et al.l 120111) . we compute the likelihood d24t 
for each SN by first introducing the hidden variables Mi, Xi, c% 
and Zi, which are, respectively, the true (unknown) values of its 
absolute B-band magnitude, stretch and colour corrections, and 
redshift. These are then assigned priors, which themselves con- 
tain further nuisance parameters, and all the introduced parame- 
ters are marginalised over to obtain the likelihood \24\ . The de- 
tails of this procedure are giv en in Appendix A of our compan- 
ion paper JMarch et al.l 120121) . By assuming separable Gaussian 
priors on the hidden variables and nuisance parameters, one can 
perform all the marginalisations analytically, except for two nui- 
sance parameters R x and R c that must be marginalised over nu- 
merically. The full likelihood function thus depends on the param- 
eters = {h, a, /?, cri n t, R x , R c }, and is 9-dimensional for the 
truncated SIS halo model and 8-dimensional for the NFW profile 
model. 



3.5 Priors on the sampled parameters 

To determine the Bayesian inference problem completely, it only 
remains to specify the prior 7r(©) on the parameters to be sam- 
pled. The choice of prior is particularly important for weak lens- 
ing analyses, since the problem is inherently underconstrained and 
therefore any prior information available is extremely useful. One 
should, however, be careful in the choice of priors not to impose 
too strong assu mptions, whic h may lead to erroneous inferences. 

Following lMarch et alJ J201 lh . we adopt the separable priors 
listed in Table [T] on each of the 'non-halo' (nuisance) parameters 
(Tint , ot, /3, R x and R c . On the halo parameters h, the priors for each 
model are listed in Table [2] and correspond to broad, conservative 
assumptions. 



SNIa population 








Dispersion of absolute magnitude 




W(- 


3, 0) on log o- int 


Lightcurve 








Stretch multiplier 


a 




U(0,4) 


Colour multiplier 


P 




W(0,4) 


Dispersion of stretch correction 


R x 


W(- 


-5, 2) on log R x 


Dispersion of colour correction 


Rc 


U( 


-5,2) onlogi? c 


Table 1. Priors on nuisance parameters, where U(a, b) denotes a uniform 


distribution between the limits a and b. 








Halo parameter 




Symbol 


Prior 


Truncated SIS halo 








Fiducial velocity dispersion (km s 


- 1 ) 




U(0, 300) 


Fiducial truncation radius (Mpc) 




r* 


U(0, 0.4) 


Exponent of cr-L power law 




V 


U(-5,5) 


Exponent of rt-cr power law 




7 


«(-5,5) 


NFW profile halo 








Concentration 




c 


U{0, 100) 


Fiducial virial radius (Mpc) 




r 200 


U{0, 0.4) 


Exponent of T200-L power law 




A 


W(-2,2) 



Table 2. Priors on halo parameters for the truncated SIS and NFW models 
respectively. 



4 SUPERNOVAE AND GALAXIES DATA SETS 

To apply the Bayesian analysis methodology described above to 
constrain the dark matter haloes of galaxies, we need observations 
both of supernovae and foreground galaxies. As in Jonsson et al. 
(2010b), we us ehigh-redshift(0 .1 < z < 1) SNIa from the 3-year 
SNLS data set dGuv et ap2010l) to constrain the properties of dark 
matter haloes of galaxies in the CFHTLS fields that intersect the 
SNIa lines-of-sight. 

The SNLS consists both of photometric and spectroscopic ob- 
servations. The photometry is obtained as part of the deep com- 
ponent of CFHTLS w ith the one square-degree imager MegaCam 
jBoulade et alj[2003h . The deep part of CFHTLS comprises four 
fields (Dl, D2, D3 and D4), each ~ 1 deg 2 in size, imaged in it*, 
g', r', i! and z filters approximately every 4-5 days during dark 
and gre y time, suitable for detecting super novae and building light 
curves dAstier et al. 2006; Guv et alj|2010l) . The spectroscopic ob- 
servations are used to determine t he nature of the supernovae candi- 
dates and measure their redshifts dHowell et alj2005l :l Brond er et al.l 
l2008l ; lBallandetal .1120091) . 



4.1 Supernovae 

Our initial sample consisted of 230 hi gh-redshift (0.1 < z < 1) 
SNIa from the 3-year SNLS d ata set fcuv et alj|20ld) . We note 
that this is 13 fewer than used bv ljonsson et alJd2010bh . since these 
SNIa wer e excluded from the SALT- 1 1 fits for various reasons de- 
scribed in IConlev et al] J201 ]]), and do not appear in the tables in 
that paper. The sample of 230 SNIa is further reduced, however, 
by edge effects and by some parts of the deep CFH TLS fields be- 
ing c overed by bright stars that have to be masked (Sulliv an et al.l 
2006). SNIa located too close to the boundary of the field or to a 
masked region are removed from the sample because of the lack of 
observat ions of foreground g alaxies. Details of this procedure are 
given in ljonsson et alJ fcOlObl) . Only 162 of the initial sample of 
230 SNIa fulfilled this selection criterion. 

To compute the predicted magnification due to the foreground 
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galaxies along the line-of-sight to each SNIa, we require an esti- 
mate of its sky location #sn and redshift 2sn- The sky location is 
obtained from the CFHTLS i'-band photometry. When available, z 
is taken to be the spectroscopic redshift of the SNIa host galaxy; 
otherwise the spectroscopic redshift of the SNIa itself is used. 

Selection cuts, S ALT-II lightcurve fitting and Malmquist cor- 
rections are made by the SNLS3 team and are already implemented 
in the supplied data files. As described in Section [3~TI the S ALT-II 
fitting algorithm is applied to each SNIa lightcurve to obtain best- 
fit estimates of: rh* B , the rest frame B-band apparent magnitude of 
the supernovae at maximumn luminosity; i\, the stretch parameter 
of the fitted light curve; and c, the colour excess in the B-band at 
maximum luminosity. These are supplemented by the covariance 
matrices C of the uncertain ties in the es t imate d lightcurve parame- 
ters, which are taken from lConlev et all d201ll) . 



4.2 Galaxies 

Our foreground galaxies are taken from the SNLS galaxy cata- 
logues in the deep CFHTLS fields. As discussed in Section [372] 
to perform our analysis we require for each galaxy: the redshift 
Zgai, sky position ga i, absolute B-band magnitude Mb, and spec- 
tral type t. The techniqu es used to obta i n thes e gal axy properties 
are des cribed in detail in ISullivan et alj d2006h and Ijonsson et al.l 
l l2010bh . 



5 APPLICATION TO SIMULATED SUPERNOVAE DATA 

In order to test our Bayesian analysis methodology, we first apply 
it to simulated SNIa data. First, SNIa photometric data in the ab- 
sence of gravitational lensing were simulated and fitted using the 
publicly available SNANA package (Kessler et al. 200 j|), in a iden- 
tical m anner to that described in our companion paper ( Mar ch et al.l 
|2012|) . In summa ry, the data were simulated to match closely the 
SNLS3 data set jGuv et al.ll201ol) by using the SNLS3 co-added 
simulation library files (which are publicly available as part of the 
SNANA package), a coherent magnitude smearing of 0.12, and 
colour smearing. The colour smearing effect, or broad-band colour 
dispersion model, implemented i n the data s i mulat ion is the EXP- 
POL model described by fig. 8 of lGuv et all d2010h. and the simu- 
lated Malmquist bias is based on fig. 14 of lPerrett et alj d2010h . 

The SNANA SNIa data simulation is a two-stage process that 
mimics the real data collection and analysis process. The first stage 
is the simulation of photometric data in accordance with the charac- 
teristic instrument and survey properties of the SNSL3 survey using 
the SNLS 3 simulation library files mentioned above. The second 
stage is the lightcurve fitting process in which the photometric data 
are fitted to S ALT-II templates to give estimates of the SNIa ab- 
solute B-band magnitude rhs, lightcurve stretch x± and colour c. 
At this lightcurve fitting stage, basic cuts are made to discard SNIa 
with a low signal-to-noise ratio and/or too few observed epochs in 
sufficient bands. After the lightcurve fitting stage we make a red- 
shift dependent magnitude correction to correct for the Malmqist 
bias; the correction is taken from a spline interpolation of table 4 in 
|Perrettetal.N2010h . 

The resulting simulated SNIa are randomly distributed on the 
sky across the four SNLS3 fields, outside of our masked regions, 
and have a redshift distribution appropriate for the SNLS3 survey. 
The final stage of our simulation process adds the lensing contribu- 
tion due to the (assumed) galaxy haloes along the line-of-sight to 
each SN. This is performed using the real SNLS galaxy catalogue, 



Halo parameter 


Symbol 


Value 


Truncated SIS halo 






Fiducial velocity dispersion 


tr* 


120 km s" 1 


Fiducial truncation radius 


f* 


0.07 Mpc 


Exponent of a-L power law 


V 


1 

:i 


Exponent of rt-o power law 


7 





NFW profile halo 






Concentration 


c 


20 


Fiducial virial radius 


r 200 


0.15 Mpc 


Exponent of T200-L power law 


A 


1 

3 



Table 3. Halo parameter values used in the generation of simulated 
gravitationally-lensed SNIa data. 



by assuming each galaxy to have a halo described by either a SIS 
or NFW profile with the parameters listed in Table|3] 

For the SIS model, the choice of the velocity dispersion cr* and 
truncation radius r„ correspond to the best-fit values obtained by 
Jonsson et al. (2010b). The value r\ = | corresponds to the Tully- 
Fisher relation for late-type galaxies, which should make up the 
majority of the SNLS sample, since nearly 95 per cent of them are 
star-forming galaxies. We also assume 7 = 0, which corresponds 
to a constant truncation radius ft = r», i.e. independent of the 
velocity dispersion and luminosity of the galaxy. 

For the NFW m odel, our choices of v alues for c and r^oo are 
motivated in part by Sprin gefet alj d2008l) . who performed a sim- 
ulation of Milky Way sized haloes and calculated the c and r2oo 
values for them. The r2oo values for these haloes are given in their 
table 1, and c value for the main sub-halo is given in their table 2. 
Typically, they find T2oo 0.2 Mpc and c ~ 10 — 16. We assume 
that most of galaxies in the SNLS sample are somewhat smaller 
than the Milky Way, and since it is well known that c is a decreas- 
ing function of mass, we assume c = 20, which yields a lensing 
signal of similar strength to that of our assumed truncated SIS halo. 
It should be noted, how ever, that the value of c in the simulations 
of Spri ngele7all (2008) varies with z and galaxy mass M, whereas 
we are assuming a fixed value for c, which is a naive assumption. 
We also assume that A = |, which corresponds to a constant virial 
mass-to-light ratio M/L oc L°. 

For an assumed set of halo properties, the lensing constribu- 
tion from the haloes along the line-of-sight to each simulated SN is 
then added to the SN apparent magnitude, and its uncertainty is left 
unchanged. 



5.1 Results for 162 simulated supernovae 

To match the real SNLS3 data to be analysed in Section [6] our 
first set of simulations contain 162 SNIa. Indeed, in this case, the 
positions on the sky of the simulated SNIa are fixed to those of the 
real SNLS3 SNIa. 

Table|4]lists the Bayesian log-evidence for each combination 
of simulation and analysis model, relative in each case to the null 
(no-lensing) model. For each simulation model, one sees that the 
analysis model preferred by the Bayesian evidence is that assuming 
no lensing. This suggests that the quantity and quality of the real 
SNLS3 data to be analysed in Section[6]are insufficent to obtain a 
detection of the lensing signal, and certainly do not allow one to 
distinguish between different halo models. 

For the four cases denoted by asterisks in Table[4] we plot the 
corresponding halo parameter constraints in Figs [TJ and [2] respec- 
tively. From Fig. [TJ one sees that, for the truncated SIS analysis 
model, the parameter constraints derived from simulated data with 
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Figure 1. ID and 2D marginalised posteriors distributions for the parameters h = {7, r;, ct», r*} of the truncated SIS halo model, derived from the analysis 
of 162 simulated SNIa data generated assuming no lensing (left) and a truncated SIS model (right). In the right-hand panel, true parameters are indicated by 
vertical lines and crosses in ID and 2D plots respectively. 




-0.5 0.5 1 5 1 0.05 0.1 0.5 1 50 1 00 - 0.1 02 



Figure 2. ID and 2D marginalised posteriors distributions for the parameters h = {A, c, ^qq} of the NFW halo model, derived from the analysis of 162 
simulated SNIa data generated assuming no lensing (left) and a truncated SIS model (right). In the right-hand panel, true parameters are indicated by vertical 
lines and crosses in ID and 2D plots respectively. 



Simulation model 

Analysis model No lensing tSIS NFW 

No lensing 0.0 ± 0.2 0.0 ±0.2 0.0 ± 0.2 

tSIS -1.8 ±0.2* -1.1 ±0.2* -1.3 ±0.2 

NFW -3.1 ±0.2* -2.6 ±0.2 -2.3 ± 0.2* 

Table 4. A log Z (log-evidence value relative to the null evidence) for the 
analysis of simulated data with 162 SNIa. Asterisks denote the cases for 
which the corresponding halo parameters constraints are plotted in Figs[T] 
and[2] respectively. 



and without a truncated SIS lensing signal are very similar, indi- 
cating that a SNLS-quality catalogue containing just 162 SNIa is 
insufficient to constrain the halo properties. In particular, we see 
that the ID marginals for <j» and r* both peak at zero (which cor- 
responds to no lensing signal). Another noteworthy feature is that, 
whereas the value of r\ is unconstrained, there is quite a strong con- 
straint restricting 7 to be negative. From dl7l l. one sees that this 
corresponds to a rather curious truncation radius scaling-law, for 
which galaxies with larger velocity dispersions have smaller trun- 
cation radii. Moreover, provided 77 is positive, this also corresponds 
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Simulation model 




Analysis model 


No lensing 


tSIS 


NFW 


No lensing 


0.0 ±0.2 


0.0 ±0.2 


0.0 ±0.2 


tSIS 


-1.7 ±0.2 


4.5 ±0.2* 


6.9 ±0.2 


NFW 


-3.1 ±0.2 


3.6 ±0.2 


7.2 ±0.2* 



Table 5. A log Z (log-evidence value relative to the null evidence) for the 
analysis of simulated data with 500 SNIa. Asterisks denote the cases for 
which the corresponding halo parameters constraints are plotted in Figs [5] 
and|4] respectively. 



to more luminous (and presumably more massive) galaxies hav- 
ing smaller truncation radii. This phenomenon was also noted by 
Jonsson et al. (2010b), but was erroneously interpreted as being a 
constraint on the halo properties derived from the SNIa data, which 
is clearly not the case since it occurs even in the analysis of simula- 
tions containing no lensing signal. Rather, from J 1 3 1 > and J 1 7b . one 
sees that a negative value of 7 allows for a smaller value of the trun- 
cation radius r t to 'offset' an increase in the velocity dispersion a, 
thereby reducing the lensing signal produced by the putative halo, 
as required to be consistent with simulated data containing no lens- 
ing signal. 

From Fig.Q] one also sees that the constraints on NFW halo 
parameters are almost identical when analysing data with and with- 
out the lensing signal. Once more, it is not possible to place con- 
traints on the halo parameters. One should note, however, that the 
exponent A of the scaling law between virial radius and luminosity 
does appear relatively well constrained around A ~ i, which corre- 
sponds to a constant mass-to-light ratio M/L oc L . Nonetheless, 
one sees a similar constraint when analysing simulated data con- 
taining no lensing signal; essentially, larger values of A (through its 
dependence on r2oo) result in a stronger lensing signal that is ruled 
out by the data. 

5.2 Results for 500 simulated supernovae 

To test that our analysis procedure is capable of detecting the grav- 
itational lensing signal and placing the correct constraints on halo 
parameters in presence of more data, we also analyse a simulated 
SNLS3-style catalogue containing 500 SNIa. In this case, the SNIa 
are randomly distributed on the sky over the four SNLS3 fields, 
outside of our masked regions. This simulated catalogue contains 
approximately three times the number of SNIa as considered pre- 
viously and is thus representative of what could be achieved by the 
SLNS programme in a total of about 9 years of observation. This 
is clearly rather unrealistic in terms of the required observing time 
and resources, but still provides a useful insight into the quantity of 
SNIa data required to make a robust detection of the lensing signal 
and constrain halo properties. 

Table|5]lists the Bayesian log-evidence for each combination 
of simulation and analysis model, relative in each case to the null 
(no-lensing) model. For the no-lensing simulation model, one sees 
that the no-lensing analysis model is again correctly preferred over 
both lensing halo models. More importantly, one now finds that for 
both simulations containing a lensing signal, the analysis models 
also containing a lensing signal are preferred at high significance, 
thereby indicating a robust detection of the gravitational lensing 
effect. Moreover, we note that in each case the correct halo model is 
preferred, although log-evidence differences between halo models 
are quite small, so this does not constitute a robust selection of the 
correct halo model. 

For the two cases denoted by asterisks in Table \5\ we plot 




Figure 3. ID and 2D marginalised posteriors distributions for the parame- 
ters h = {7, r), cr» , r* } of the truncated SIS halo model, derived from the 
analysis of 500 simulated SNIa data generated assuming a truncated SIS 
model. True parameters are indicated by vertical lines and crosses in ID 
and 2D plots respectively. 

the corresponding halo parameter constraints in Figs [3] and [4] re- 
spectively. From Fig. [3] which corresponds to a truncated SIS halo 
model in both the simulation and the analysis, one sees that each 
of the ID and 2D marginalised posteriors has a well-defined peak 
(away from zero), indicating a robust detection of the lensing sig- 
nal. This supports the conclusion drawn above from the Bayesian 
evidence results in Table [5] One also sees that all the parameters 
have been estimated with reasonable accuracy, although some of 
the 2D marginals exhibit pronounced degeneracies between the pa- 
rameters, particular those involving the parameter r„, which is the 
least well constrained. 

Similar conclusions can be drawn from the parameter con- 
straints corresponding to a NFW halo model in both the simulation 
and analysis, which are shown in Fig|4] One sees that the ID and 
2D marginalised distributions have reasonably well-defined peaks 
that contain the true parameter values, although again degeneracies 
do exist between parameters. One also sees that the ID posteri- 
ors for both c and r^oo peak well away from zero, supporting the 
conclusion drawn from the Bayesian evidences in Table|5]that one 
obtains a strong detection of the lensing signal. 



6 APPLICATION TO REAL SUPERNOVAE DATA 

We now apply our Bayesian analysis methodology to the real SNIa 
and galaxies data sets described in Section [4] For the real data 
we again consider both the truncated SIS and NFW halo models 
in our analysis. For each case, we also investigate dividing the 
foreground galaxies in the SNLS catalogue into passive and star- 
forming classes, each having halo model parameters that are al- 
lowed to be independent. 

Table [4] lists the Bayesian log-evidence for each analysis 
model, relative in each case to the null (no-lensing) model. In the 
case of one galaxy type, one sees that there is a very marginal pref- 
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Figure 4. ID and 2D marginalised posteriors distributions for the param- Figure 5. ID and 2D marginalised posteriors distributions for the parame- 

eters h = {A, c, r| 00 } of the NFW halo model, derived from the anal- ters h = {7, 77, <r» , r„} of the truncated SIS halo model, derived from the 

ysis of 500 simulated SNIa data generated assuming a NFW model. True analysis of real SNLS3 data, 

parameters are indicated by vertical lines and crosses in ID and 2D plots 
respectively. 



1 galaxy type 2 galaxy types 

Model (faP = h sC ) (h p £ h sl ) 

Truncated SIS 0.2 ±0.2* 0.5 ±0.2* 

NFW -2.5 ±0.2* -2.2 ±0.2 

Table 6. A log Z (log-evidence value relative to the null evidence for no 
lensing signal) for the analysis of the real SNIa data. h p = h s{ indi- 
cates that the halo parameters for passive and star-forming galaxies are 
constrained to be the same, whereas for h p 7^ h si they are allowed to be 
independent. Asterisks denote the cases for which the corresponding halo 
parameters consraints are plotted in Figs|3]|6]and|7] respectively. 

t 

erence of just 0.2 log-evidence units for a lensing signal produced 
by a truncated SIS halo model. This is, however, only at the level 
of the uncertainty in the evidence calculation, and so this model 
and the no-lensing model are broadly equally favoured. A lensing 
signal produced by a NFW halo model is, however, strongly dis- 
favoured with an evidence of —2.5 log-evidence units relative to 
the no lensing model. 

When one splits the foreground galaxies into their passive and 
star-forming spectral types, one sees that the evidence for a lensing 
signal produced by a truncated SIS increases very slightly and is 
marginal preferred over the no lensing model by 0.5 log-evidence 
units. Once again, a lensing signal produced by a NFW halo model 
is quite strongly disfavoured by —2.2 log-evidence units. 

We discuss the halo parameters constraints for these various 
models in the following three subsections. 

6.1 Truncated SIS model 

The ID and 2D marginalised parameter constraints obtained for the 
truncated SIS halo model with one galaxy type are shown in Fig.[5] 
As one might expect from the relative evidence for this model, the 



parameters constraints lie mid-way between those obtained for the 
simulated SNIa data with 162 and 500 supernovae respectively, as 
discussed in Section[5] 

Focussing on the ID marginal posterior distribution of each 
halo parameter, we first note the strong preference for negative val- 
ues of 7, as we observed in the analysis of simulations containing 
no lensing signal. As explained in Section|5] this is merely a conse- 
quence of the lensing signal (if any) in the data being very weak and 
does not constitute a meaningful constraint on the halo properties. 

Conversely, there is a relatively strong preference for positive 
values of rj, corresponding to a positive correlation between the 
halo velocity dispersion and luminosity. Indeed, the marginal dis- 
tribution peaks at rj ~ 1, which is considerably larger than the 
canonical values of rf — i and r] si — |, corresponding to the 
Faber-Jackson and Tully-Fisher relations, valid for early and late- 
type galaxies, respectively. This constraint also differs noticeably 
from that shown in the left-hand panel of Fig. Q] obtained from 
simulated data containing no lensing signal, which is suggestive of 
a marginal lensing signal being present in the real data. 

The possibility of a faint lensing signal being present in the 
data is also suggested by the ID marginal distribution for the fidu- 
cial velocity dispersion <r* , which has a clear peak away from zero. 
Indeed, the peak is centred on <r» ~ 120 km s , which was also the 
value assumed in the simulations analysed in Section [5] The fact 
that the parameter constraint obtained from the real data again ap- 
pears tighter than that obtained from the simulations of 162 SNIa 
with no lensing signal, shown in the left-hand panel of Fig.Q] may 
be a result of the real data preferring 77 ~ 1, as opposed the value 
rj = 1 used in our simulations. From J 1 3 1 > and d!7l >, one sees that 
the lensing signal produced by a halo with a given cr* is larger in 
the former case. 

Finally, we note that the ID marginal for the fiducial trunca- 
tion radius peaks very close to zero and yields no real constraint on 
r*, as we found in our analysis of the 162 simulated SNIa. 
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Figure 6. ID and 2D marginalised posteriors distributions for the parame- 
ters h = {A, c, ?*20o} °f me NFW halo model, derived from the analysis of 
real SNLS3 data. 

6.2 NFW model 

The ID and 2D marginalised parameter constraints obtained for the 
NFW halo model for one galaxy type are shown in Fig. [6] The con- 
straints on the parameters are very similar to those shown in Fig[2] 
obtained from 162 simulated SNIa with and without a lensing sig- 
nal produced by a NFW halo model with the parameters listed in 
Table [3] One sees that it is not possible to place meaningful con- 
straints on any of the parameters, indicating the absence of a de- 
tectable lensing signal. This is consistent with the low relative evi- 
dence for this model listed in Table|6] 

6.3 Separating passive and star-forming galaxies 

So far we have assumed that the halo parameters for the trun- 
cated SIS and NFW models, respectively, are the same for differ- 
ent galaxy types. This may be an oversimplication and so it is of 
interest to split the foreground galaxies into their passive and star- 
forming spectral types, both of which are allowed to have indepen- 
dent halo parameters. 

We concentrate on the truncated SIS model, since this yielded 
a significantly higher evidence than the NFW profile in our analysis 
assuming just a single galaxy type (see Table[6]l. The galaxies in our 
sample are classified as either passive or star-forming depending on 
their sSFR (see Section POt . We thus allow all the halo parameters 
to be different for each type of galaxy, so that the full halo param- 
eter space becomes h = {7 P , rf , uf , rj, 7 sf , ?7 sf , , rj f }. The 
fraction of passive galaxies ranges from 5 per cent in the CFHTLS 
fields Dl and D3 to 8 per cent in field D4. Thus, the vast majority 
of foreground galaxies are star-forming. 

The resulting ID and 2D marginalised posteriors for the pa- 
rameters are shown in Fig. [7] One sees that the constraints on the 
parameters for passive and star-forming galaxies are very differ- 
ent. As one might expect, given the relative percentages of passive 
and star-forming galaxies in the catalogue, the constraints on the 
halo parameters of the passive galaxies are very much weaker than 



those for star-forming galaxies. Indeed, for the passive galaxies, 
the marginal distributions closely resemble those shown in Fig.fj] 
obtained from simulations containing 162 SNIa with no lensing 
signal. Conversely, for star-forming galaxies, the constraints are 
slightly tighter than those shown in Fig. [5] obtained from the real 
data assuming just one galaxy type, but resemble them in their main 
features. In particular, we again see a strong preference for negative 
values of 7 s , but now the ID marginal distribution has a modest 
peak at 7 sf ~ —3. We also recover a strong constraint that rj af is 
positive, with the ID marginal possessing a small peak at r/ si ~ 1. 
The ID marginal for af peaks strongly away from zero, preferring 
a value of af ~ 150 km s _1 . There is, however, no constraint on 
rj , which is in keeping with our previous findings. 

Finally, it is worth noting that the 2D marginal in the 
(trf , al )-space shown in Fig.|7jb ears a passing resemb lance to the 
corresponding plot (figure 6) in I Jonsson et al] i feOlObh , but only if 
a, and af are interchanged. Given the relative percentages of pas- 
sive and star-forming galaxies in the catalogue, it seems sensible 
that one should obtain a tighter constraint on af than cr», as we 
find in Fig. [7] By contrast, figure 6 in Jonsson et al. (2010b) shows 
the opposite, which suggests that they have erroneously swapped 
the galaxy types in their analysis. 

6.4 Inclusion of additional low-redshift supernovae data 

We also investigate the inclusion of additional, low-redshift SNIa 
data to our analysis as a potential means of enhancing the detec- 
tion of a halo lensing signal. The rationale here is first to anal- 
yse the low-z SNIa data alone, using precisely the same method- 
ology as for the SNLS3 data, but assuming no lensing from fore- 
ground galaxy haloes. The resulting posterior distributions derived 
for (most of) the nuisance parameters listed in Table Q] may then 
used as priors on these parameters in the subsequent analysis of the 
(high-z) SNLS3 data, thereby replacing the very conservative pri- 
ors assumed in TableQ] which may be 'diluting' the lensing signal. 

As our low- z SNIa dataset, we use an updated version of the 
catalogue compiled by Sullivan et al. (2010). This consists of a 
sample from the compilation of Conley et al. (2010), which itself 
includes SNIa from a variety of sources (primarily Hamuy et al. 
1996a; Riess et al. 1999; Jha et al. 2006; Hicken et al. 2009a; Con- 
treras et al. 2009). Sullivan et al. apply bulk-flow peculiar velocity 
corrections to the SN magnitudes and redshifts, placing the red- 
shifts in the CMB-frame (z cmb ) following Neill et al. (2007), but 
with updated models (Conley et al. 2010). Only SN la in the smooth 
Hubble flow, defined as z cm b > 0.01 are used, and the same light 
curve quality cuts as for the SNLS3 sample are employed. There 
are 123 low-redshift SN la in total, with redshifts in the range 
[0.01,0.08]. 

Analysis of these data yields approximately Gaussian pos- 
teriors (truncated to positive values) on the nuisance parameters, 
given by: a ~ A/"(1.30, 0.11 2 ), /3 ~ Af(2.81, 0.18 2 ), R x ~ 
A/"(0.119,0.008 2 ) and R c ~ Af(0.Q74, 0.05 2 ). These distribu- 
tions are clearly much tighter than the original assumed priors given 
in Table Q] note that we use the original prior on the nuisance 
parameter <Ti nt , as this quantity is likely to differ between low-2 
and high- z SNIa, since observational uncertainties such as the es- 
timation of photometric errors can lead to variations of amt on a 
sample-by-sample basis dMarch et alj2012b . We find that using the 
posteriors derived from the low-z SNIa as priors in the subsequent 
analysis of the high-z SNLS3 data yields parameter constraints that 
are almost identical to those presented in Figs [5j[7] Moreover, the 
log-evidences relative to the null (no-lensing) model are found to 
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be very similar to those listed in Table|6] This finding makes sense, 
since the new, tighter priors on the nuisance parameters are used 
for both the halo and null (no-lensing) models, and they are, in fact, 
consistent with the posteriors derived on these parameters from the 
SNLS3 data alone. Thus, in summary, our results are unchanged by 
the inclusion of low-z SNIa data. 



7 CONCLUSIONS 

We have presented a Bayesian statistical methodology for con- 
straining the properties of dark matter haloes of foreground 
galaxies that intersect the lines-of-sight towards SNIa. The 
method builds upo n the Bayesian hierarchical model presented by 
iMarch et al.l d201 ll) for improving constraints on cosmological pa- 
rameters from SNIa observations. Compared with the usual \ 2 ' 
method, which suffers from shortcomings in term s of its statistical 
foundations and robustness, IMarch et alj J201 lh demonstrate that 



the Bayesian method delivers tighter statistical constraints, reduces 
statistical bias and produces confidence intervals with better statis- 
tical coverage. 

We use this methodology to obtain an effective likelihood 
function giving the probability of obtaining the observed SNIa data 
(i.e. the parameter values obtained in SALT-II lightcurve fits) as 
a function of the parameters of the dark matter halo models as- 
sumed for the galaxies along the lines-of-sight to the SNIa. Follow- 
ing the imposition of suitable priors on these parameters (together 
with some nuisance parameters), we explore the full posterior dis- 
tribution in all the parameters simultaneously using the nested sam- 
pling algortihm MULTlNEST, which also calculates the Bayesian 
evidence for use in model comparison. 

We first apply our method to simulated SNIa datasets gen- 
erated using 162 high-redshift (0.1 < z < 1) SNIa from the 
3-year data release of the Supernova Legacy Survey (SNLS3) as 
a template and assuming either a truncated singular isothermal 
sphere (SIS) model or a Navarro-Frenk-White (NFW) model for 
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the dark matter halo density profile of foreground galaxies in the 
deep Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) 
fields. These simulations were generated using the SNANA pack- 
age, assuming realistic values for the halo parameters and obser- 
vational data quality. Assuming conservative priors on the parame- 
ters, we demonstrate that one cannot obtain a detection of the lens- 
ing signal in this case. By considering instead a simulation cat- 
alogue containing 500 SNIa, we show that our method produces 
posterior distributions for the parameters that are fully consistent 
with their input values in the simulations. Moreover, in each case, 
the Bayesian evidence for the correct class of halo model far ex- 
ceeds that of the no-lensing model, indicating a strong detection of 
the simulated lensing signal. We also analyse each simulation as- 
suming the other (incorrect) class of halo model and show that the 
Bayesian evidence is able to select the correct model in each case, 
albeit only marginally. 

In the analysis of real SNLS3 data (consisting of 162 SNIa) 
we find, contrary to previous studies, only a very marginal detec- 
tion of a lensing signal in the case of the truncated SIS halo model, 
which is preferred by just 0.2 log-evidence units relative to the no 
lensing model. Indeed, since this difference is similar to the uncer- 
tainty in the evaluation of the evidence, one may consider the no 
lensing model to be equally favoured by the data. Nonetheless, the 
parameter constraints for the truncated SIS halo model do appear 
somewhat tighter than those obtained for simulations of 162 SNIa 
without the inclusion of a lensing signal, which is again sugges- 
tive of a borderline detection of a lensing signal in the real SNLS3 
data. One finds that the SNLS3 data strongly prefer negative val- 
ues of 7, which corresponds to luminous galaxies having smaller 
truncation radii than less luminous ones, but this is simply a mani- 
festation of the lensing signal being very weak and does not consti- 
tute a meaningful constraint on the halo properties. There is also a 
preference for positive values of rj, corresponding to a positive cor- 
relation between the halo velocity dispersion and luminosity, with 
the ID marginal for this parameter peaking at r\ « 1, which is 
somewhat larger than the canonical values of rf — j and r/ sl = | 
valid for early and late-type galaxies, respectively, and leads to a 
stronger lensing signal for a halo with a given luminosity. The ID 
marginal distribution for the fiducial velocity dispersion <t» has a 
clear peak away from zero, which is centred on <r» ~ 120 km s , 
but the marginal for the fiducial truncation radius peaks very close 
to zero and yields no real constraint on r, . No detection or mean- 
ingful halo parameter constraints are obtained for the NFW model. 

Finally, we investigate the possibility that the halo parameters 
may be different for passive and star-formation galaxies, respec- 
tively. Focussing on the truncated SIS model, we find that allowing 
all the halo parameters for the two galaxy types to be completely 
independent increases the evidence for the model slightly, by 0.3 
log-evidence units. Since only around 5 per cent of the galaxies in 
our catalogue are passive, we find that one cannot place any mean- 
ingful constraints on their halo properties, but that the constraints 
on the halo parameters of the star-forming galaxies are similar to 
those obtained if one assumes just one galaxy type. 

Our results contradict to some extent the previously re- 
ported high-significance detections of gravitational lensing of SN Ia 
(Jonss on et al J2007I : iKronborg etaifeoid : I Jonsson et all2010al fbb. 
where the last study uses essentially the same SNLS3 data as those 
analysed here. We have also verified that our findings are un- 
changed by the inclusion of additional low-z SNIa data. The major 
difference between these earlier analyses and the study presented 
here is the s t atisti cal methodology employed. As demonstrated by 
March ^tal] ( l201ll) . the usual x 2 -method used in previous analy- 



ses has a number of shortcomings in terms of its statistical foun- 
dations and robustness, including not allowing for rigorous model 
checking and not providing a reliable framework for the evalua- 
tion of systematic uncertainties. This should be contrasted with the 
statistically-principled and rigorous Bayesian hierarchical model 
used here. Moreover, previous studies employed using simple grid- 
ding methods to explore the parameter space of the dark matter 
halo models, which allow only a small subset of the parameters to 
be varied simultaneously. 

It should be mentioned, however, that the analysis presented 
here does make the significant simplifying assumption that all dark 
matter haloes of the foreground galaxies have the same values for 
their free parameters (at least within their spectral type). This is not 
too restrictive an assumption for free parameters corresponding to 
fiducial values appearing in scaling laws, but it is certainly a sub- 
stantial oversimplication for physical parameters such as the con- 
centration c in the NFW profile. This may weaken the constraints in 
the analysis of real data, as compared with the analysis of simulated 
data generated assuming all the haloes to have the same parameter 
values. We will investigate this possibility in a future work 

In closing, we note that the analysis code used in this paper is 
quite general in nature and could easily be applied to other combi- 
nations of SNIa and foreground galaxy catalogues. Anyone wishing 
to use the code in collaboration should contact the authors. 
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